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Preliminary results of identical-particle correlations probing the geometric substructure 
of the particle-emitting source at RHIC are presented. An mT-independent scaling of 
pion HBT radii from large (central Au+Au) to small (p+p) collision systems naively 
suggests comparable flow strength in all of them. Multidimensional correlation functions 
are studied in detail using a spherical decomposition method. In the light systems, the 
presence of significant long-range non-femtoscopic correlations complicates the extraction 
of HBT radii. 

1. Introduction 

Particle interferometry is a useful technique that provides information on the space- 
time properties of nuclear matter created in high energy collisions (for the latest review 
articles see [II]I2])- In this paper results of particle correlations in p-|-p, d-l-Au and Cu-|-Cu 
collisions at ^/sNN = 200 GeV and Au+Au at ^/sNN = 200 and 62 GeV registered by 
the STAR experiment are presented. Rich data statistics and the wide acceptance of the 
detectors gives an opportunity to do a multidimensional femtoscopic analysis. The main 
focus of this article is on the transverse mass (m^ = + ) dependence of the HBT 
radii for different system sizes. This allows the study of the dynamics of the source and the 
conditions of the collision. Importantly, direct comparison between the mj- dependence 
of radii from p+p and A+A collisions is possible for the first time. 

2. Multiplicity Scaling and Expansion 

Figure m presents the HBT radii dependence on (dNch/dr]Y^^ {dNch - number of charged 
particles) for different colliding systems at different energies of the collisions. The motiva- 
tion for studying such a relation is its connection to the final state geometry through the 
particle density at freeze-out. All STAR results, from p+p, d+Au, Cu+Cu and Au+Au 
collisions, are combined on the left panel of this figure and, as seen, all radii exhibit a 
scaling with {dNch/dr]Y^^. On the right panel STAR radii, this time for different range 
of kx, are plotted together with AGS/SPS/RHIC systematics [ IJ. It is impressive that 
the radius parameters Rside and Riong follow the same curve for different collisions over 
a wide range of energies and, as it was checked, this observation is valid for all /ct bins 
studied by STAR. It is a clear signature that the multiplicity is a scaling variable that 

*For the full list of STAR authors and acknowledgments, see appendix 'Collaborations' of this volume. 
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Figure 1. Femtoscopic radii dependence on the 
number of charged particle. 
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Figure 2. Final size of the source vs. ini- 
tial radii calculated from Glauber model. 
STAR Au+Au data from [gl. 



drives these HBT radius parameters. Rout mixes space and time information. Therefore 
it is unclear whether to expect a simple scaling with the final state geometry. 

Despite the approximately linear dependence of Rside and Riong on d{dNch/ drjY^^ , be- 
cause of the finite intercepts of the linear relationships [d EI, the data do not imply 
freeze-out at a constant density. The fact that the radii scale at all (linearly or not) 
indicates that the sideward and longitudinal radius is determined only by multiplicity, 
independent of energy, colliding system, or impact parameter [IHEIC]. This scaling, how- 
ever, breaks down at much lower energies, when baryons constitute a significant fraction 
of the freeze-out system [ IHl C] . 

Figure 121 suggests that the relationship between initial and final geometry is not trivial, 
however. There, the final RMS of the source, estimated by Rside at low /cr [Ej, is plotted 
versus the initial RMS of the overlap region, estimated with a Glauber calculation. The 
system generated in central Au+Au collisions undergoes a two-fold expansion, while those 
from the most peripheral d+Au and p+p collisions expand little. Does this imply that 
small systems are less explosive than large ones? This is explored in the next Section. 



3. Transverse Mass Dependence of HBT Radii : p+p vs Au+Au 

The iriT dependence of the femtoscopic radii in heavy ion collisions is usually attributed 
to collective flow of a bulk system, and recent STAR results for central Au+Au collisions 
support this picture [|H|. In a flow scenario, an approximately "universal" m-r dependence 
should apply to all particle types. This is in fact observed in Figure El in which one- 
dimensional radii from pion [|^, charged kaon [in|, neutral kaon proton and anti- 
proton [E], and proton- A [ 12j correlations are plotted. Correlations between particles 
with very different masses also show characteristic signals of collective flow [IT^. 
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Figure 3. m-p dependence of Rinv for dif- 
ferent particles. PHENIX data from [ 9 . 



Figure 4. Ratio of HBT radii from Au+Au, 
Cu+Cu and d+Au by p+p collisions at 
y/sj^ = 200 GeV. 



Decreasing ttit dependences of the radii have previously been reported for hadron [ IT^ 
and e"*" + e~ [ |^ collisions. Until the STAR analysis, a direct comparison of the rriT 
systematic from hadron collisions with those from A+A collisions had been hampered by 
the use of quite different parameterizations in the different experiments. As shown in 
Figure |3J the ratio of HBT radii from Au+Au collisions with those from p+p collisions 
are approximately independent of my. It might be that proposed explanations for the mr 
dependence in small systems, including Heisenberg-type relations [^j and inside-outside 
cascade scenarios [ CSj, replicate the rriT dependence from bulk collective flow. This 
would be a surprising and unfortunate coincidence, suggesting that dynamical signals in 
the space-time sector cannot distinguish very different underlying physics. Alternatively, 
as proposed by Csorgo et al. [ 17J, it could be that p+p collisions generate a thermalized 
hulk system, similar to that created in heavy ion collisions. Further theoretical work on 
this important question would be most welcome. However, as it is mentioned in the next 
Section, further work is also required on the experimental side. 



4. The effect of non-femtoscopic correlations in low-multiplicity events. 

HBT radii may arise from Gaussian fits to the three-dimensional correlation function [ 
p. But for low- multiplicity systems, unaccounted-for structures in the correlation func- 
tion itself are observed. In particular, the value of the correlation function C{q) for large 
\q\ (larger than the scale of quantum statistics or Coulomb interactions) does not approach 
a common g^independent value. Thus, non-femtoscopic correlations (e.g. due to momen- 
tum conservation) are becoming significant here; they have been also observed earlier, in 
elementary particle collisions, see e.g. [I14j. 

The spherical decomposition of the correlation function in g^space has been proposed [ 
ITS] as a sensitive measure of long-range non-femtoscopic correlations. The lowest-/ com- 
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ponents Ai^m of this procedure are shown in Figure El For all correlations of femtoscopic 
origin (Gaussian or not), A;^o,m must vanish [HH] for large q; clearly other correlations 
are present in STAR data. 

While various possible sources of these corre- 
lations are still being explored, it is interesting to 
simply account for the long-range correlations by 
adding ad-hoc terms to the Gaussian fits, which 
produce constant values of A2fi and ^42,2- It was 
found that the ratios shown in Figure 4 remain 
mT-independent, though their values shift some- 
what. 
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5. Conclusions 
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Figure 5: First three non- vanished 
components of the spherical harmonic 
decomposition of the correlation func- 
tion for peripheral d+Au collisions. 



The results of pion interferometry for all en- 
ergies and colliding systems at RHIC have been 
presented. In agreement with data at SPS and 
AGS, STAR indicates that the multiplicity is the 
scaling variable that determines the size of the 

source at freeze-out. Perhaps surprisingly, the rriT dependence of HBT radii appears to 
be independent of collision species or multiplicity. Finally, a problem with the baseline of 
the correlation function for low multiplicity collisions has been reported, and a promising 
tool based on the spherical harmonic decomposition of the correlation function has been 
used in order to address it. The physics of this structure remains under investigation. 
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